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Abstract 

The role of graphene in mediating the exchange interaction is theoretically investigated when it 
is placed between two ferromagnetic dielectric materials. The calculation based on a tight-binding 
model illustrates that the magnetic interactions at the interfaces affect not only the graphene band 
structure but also the thermodynamic potential of the system. This induces an indirect exchange 
interaction between the magnetic layers that can be considered in term of an effective exchange 
bias field. The analysis clearly indicates a strong dependence of the effective exchange bias on 
the properties of the mediating layer, revealing an effective mechanism of electrical control even 
at room temperature. This dependence also results in qualitatively different characteristics for the 
cases involving mono- and bilayer graphene. 

PACS numbers: 75.70.Cn,75.80.+q,75.75.+a,73.21.Ac 
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The exchange bias is a phenomenon associated with the unidirectional anisotropy of the 
exchange interaction at the interface of ferromagnetic (FM) and antiferromagnetic (AFM) 
materials Widely used in the magnetic sensors, magnetic random access memory, read 
heads of hard disk drives, etc., its most important characteristic is the exchange bias field 
that shifts the centrum of the hysteresis loop. Observation of this phenomenon has generally 
been attributed to the use of compounds whose FM component has a Curie temperature 
higher than the Neel temperature of the AFM counterpart. Hence the structure, when cooled 
in the presence of an external magnetic field, first aligns the FM layer which in turn orders 
the AFM component accordingly. In the proximity of an AFM pinning layer at or below the 
Neel temperature, the FM material thus experiences the unidirectional anisotropy induced 
by the exchange interaction. The temperature dependence can also lead to the unusual effect 
of spontaneous magnetization reversal through an interplay between the external magnetic 
field and the exchange bias field aligned in the antiparallel direction ^|. In some cases, the 
exchange bias field may degrade gradually after multiple recycling of the system through 
the consequential hysteresis loops (i.e., the so-called training effect in the polycrystalline 
magnetic materials) []], [3, Q] . 

As all the properties described above are related to the magnetic phenomena or the 
magnetic field repercussion, it raises a natural question: Can an electric field control the 
exchange bias field at the magnetic interfaces? One potential solution may be to utilize 
the multiferroic films in place of the FM layers as it has been discussed recently in Ref. jsj. 
However, the varietv of multiferroic materials that can operate at room temperature have so 
far been limited [sl,lf3]. On the other hand, a completely different approach may be possible 
by exploiting the unique properties of atomically thin graphite (i.e., graphene) when placed 
at the interface between two magnetic dielectric layers. Previous analyses demonstrated 
that the novel phenomena induced by the exchange interaction between FM layers and 
monolayer graphene (MLG) 7, 8] or bilayer graphene (BLG) depend sensitively on the 
detailed electronic features. It certainly offers a ground to explore electrical control of the 
magnetic interactions, for which the graphene electrons play a vital role. 

In this work, we explore the exchange interaction of dielectric magnetic layers mediated 
by graphene and its manipulation by electrical control. Our theoretical result clearly illus- 
trates that this interaction can be expressed in terms of the effective exchange bias field. It 
also reveals that the strength of the effective bias field depends critically on the electronic 
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properties of graphene, particularly the position of the electro-chemical potential \i (i.e., the 
carrier density), amenable for electrical control [h]]. Numerical estimates of the predicted 
phenomenon are provided for an order of magnitude analysis, illustrating its potential sig- 
nificance in spintronic applications. 

For accurate understanding, it is useful to identify the underlying differences between 
those involving MLG and BLG. In the case of MLG, the electron spins polarized by one FM 
layer interact with the other FM layer mediating their parallel alignment. Therewith the 
indirect interaction between FM layers amplifies with an increase in the electron density in 
MLG. On the other hand, the two carbon layers in BLG interact primarily with the nearest 
magnetic dielectric, respectively, leading to the possibility of inter-layer symmetry break 
through nonidentical exchange interactions [9|. For example, the difference in the interaction 
can be achieved by varying the mutual orientation of FM layers that are in contact with 
BLG at the top and bottom interfaces. In addition to the modified BLG band structure as 
discussed in Ref. [9[ (similar to the case of an external electrical bias III, ll2, HQ 
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the impact of this asymmetry is .note fundamental by energetieally favoring an antipatallel 
configuration in the magnetization orientation of two magnetic layers [16( . Hence, both MLG 
and BLG can be considered as mediating the exchange interaction (and, thus, the exchange 
bias) between FM layers, each with a different preference for the bias field orientation. 

We begin with a brief analysis of the graphene electronic band structure sandwiched be- 
tween two magnetic dielectrics. The magnetization belonging to each dielectrics, or more 
precisely to the proximate strata of FM (or AFM) films in contact with graphene is 
defined as Mi and M2, respectively. A particularly interesting feature appears for ferro- 
magnetic dielectric layers (FDLs) with an identical magnitude of the magnetic moments 
(| Mi I = IM2I = M) so that the total magnetization can be varied from zero to 2M. If 
the absence of electrically induced asymmetry is further assumed, then only the relative 
misalignment of Mi and M2 provides the cause of symmetry break between two graphene 
layers in the case of BLG. 

Figure 1 schematically illustrates the specific structure under consideration. It resembles 
the ferromagnet-metal-ferromagnet hybrid structures that reveal a giant magnetoresistance 
owing to the spin-dependent in-plane conductivity [181 ] . In our case, graphene substitutes 
the metallic film. The magnetization M x of the bottom FDL is assumed to pinned along 
the x direction by an AFM substrate with a sufficiently strong exchange bias field. The top 
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FDL is constructed from the same material but its magnetization vector M2 can be rotated 
in the x-y plane (by an external or "internal" magnetic field) forming an angle 9 with Mi. 

The influence of the FDL magnetization on the graphene electronic structure can be 
realized in actual structures through either the exchange interaction with magnetic ions 
(assuming an overlap between the carbon 7r-orbitals and unfilled shells of the magnetic ions 
in the FDLs) or an interaction via the ligands of FDLs. Thus, the problem can be modeled 
in the mean field approximation 7|, |8| with the Hamiltonian 

where TIq is the spin-independent MLG (n = 1) or BLG (n = 2) Hamiltonian. The 
remaining term of Eq. ([1]) describes the energy of an electron spin S in the effective fields (in 
units of energy) aMj and qM 2 of the proximate FDLs, where parameter a is proportional 
to the carrier- ion exchange constants as evaluated in Refs. [zj and sj]. Hence, it can be 
written for MLG as 

W« =a(Mi + M 2 )S, (2) 



while the corresponding expression for BLG is 

ftg) = V 1 aM 1 S+V 2 aM 2 S . (3) 

In Eq. (J3J), the projection operator V\ (V2) is 1 for the electron localized at the bottom (top) 
carbon monolayer and otherwise. 

In the case of low energy electronic excitations, the tight-binding approximation accu- 
rately describes the band spectra of graphene in the vicinities of each valley K and K' 



12 . Il9j |. Hence, we adopt the tight-binding Hamiltonian near the valley extrema aug- 
mented to account for the electron exchange energy in the form of Eq. (jSJ) or Eq. Q. The 
qualitative difference between the two spin Hamiltonians [i.e., Eqs. (j2J) and d3J)] implies 
that Tt^J = 2aM cos(9/2)Sm (Sm is the electron spin projection on the direction of sum 
Mi + M2) commutes with h}q while this is not the case for H^x and TCq . As a result, 
MLG remains gapless with four non-degenerate branches e^(k) near the K and K' points 
(except the case of 6 = rr that corresponds to Tiex = 0). They are identical for the conduc- 
tion (b = +1) and valence (b = —1) bands (i.e., the reflection symmetry about e = 0) and 
isotropic with respect to the valley centrum k — 0. The degeneracy of spin doublet is lifted 
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as stated, which leads to the expression 

Q 

£ W (k) = bhv F k + oaM cos - . (4) 

Here, vp = 10 8 cm/s is the Fermi velocity in MLG, the subband index a = ±1, and 9 is the 
angle between Mx and M 2 (0 < 9 < n). The impact of the spin Hamiltonian on the MLG 
band structure is clearly illustrated in Figs. 2(a) and 2(b) for two opposite cases of M 1( /M 2 
alignment. An important point to note is that, while MLG remains gapless independent of 9 
and aM following the qualitative discussion given above, the contact (i.e., zero gap) points 
are shifted to the circle of radius k c = (\a\M/hvF) cos#/2 in the k x -k y plane (measured from 
the centrum of the K or K' point). Once 9 = n and Mi + M2 becomes zero in Eq. (T5]), the 
band structure returns to that of unaltered MLG as shown in Fig. 2(b). Consequently, the 
density of states for MLG interacting with FM layers modifies to 

Ji) fe \ _ o max{|e|,HMcosf} 

that is non-zero even at e = 0. 

For BLG with the Hamiltonian TiS 2 \ a qualitatively different situation is realized due 
to the interlayer electron transitions with the corresponding matrix element 71 = 0.4 eV. 

Following the approach discussed in Ref. [9|, the BLG energy spectra can be obtained in 

(2) 

terms of eight energy branches e b x CT (k) for each valley. Along with b and a described above, 
an additional index A = ±1 is introduced to distinguish four low-energy bands e b _ la (k) 

> 71. As it is convenient to 
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from the other four excited states with energies 
normalize the parameters in units of 71, the dimensionless momentum p = hvpk/'yi and the 
effective field G = aM/71 are used hereinafter. Then, the energy bands can be expressed 
as 

4 2 !» = hi^P 2 + ^ + \ (l + ^Gcos 6 - + XWM^j , (6) 

where 



9\ 2 , „ 4 / „ 9^ 2 



W a (p) = \l [1 + aGcos-j (l + 4p 2 )+ (^pGsin-J . (7) 

Four solutions with £+i,a,ct > correspond to the conduction bands, while their mirror 
images with respect to the zero energy describe the highest valence bands (b = —1). 

In contrast to MLG, the calculation clearly illustrates the presence of an energy gap 
Eg between the lowest conduction band and the highest valence band as soon as 9 ^ 0. 



When the orientation of Mi and M2 is in a parallel alignment (9 = 0), the net effect of 
the exchange interaction simply lifts the two-fold spin degeneracy resulting in two pairs 
of spin-split bands that cross each other at p = a/ G(l + G/2) /2 as shown in Fig. 2(c). 
However, once they are misaligned, the electronic bands become of mixed spin character 
(e.g., with both parallel and antiparallel components to the x direction). Subsequent anti- 
crossing opens up the gap that progressively grows with 9. At 9 = it in Fig. 2(d) (i.e., 
Nli — — M2), the gap reaches the maximum while the energy bands regain the two-fold 
degeneracy similar to MLG [9]. Accordingly, the BLG density of states p^\e) is zero for 



kl < £9/2, where e g = 71 G sin §/y 1 + G 2 + 2G cos |. As the density of states determines 
the characteristic carrier occupancy, the calculation result strongly indicates that the total 
electronic energy of the structure can be controlled by the orientation of Mi and M2 or, 
more specifically, by the angle 9. This point is schematically illustrated in Fig. 2 through 
the comparison of shaded regions representing the occupied states when the electro-chemical 
potential \i (dashed line) is hypothetically located at e — 0. 

Now we can take into account the entropy effects at finite temperature along with the 
band modification [Eqs. (jlj) and (jHJ)] by evaluating the thermodynamic potential 



n n (9) = -2k B Tj2J2 ln 

{m} k 



exp 



k B T 
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where k B is the Boltzmann constant, T is the temperature, n = 1, 2 denotes the cases of 
MLG and BLG, respectively, and {m} collectively represents the band indices described 
earlier. Additionally, the factor of 2 comes from the summation over the valleys K and K' . 
One remarkable outcome from the calculations of thermodynamic potential [Eq. (JH])] is the 
universal scalability of the dependence on angle 9 for both MLG and BLG: 

&n n (9) 1 , 

where AQ n (9) = Q n {9) — Q n (0). In other words, the relative change in the thermodynamic 
potential [Eq. ([9])] is independent of all parameters except cos 9 in spite of its seemingly 
complex expression [see Eqs. (jl]), and ([HD]- 

Clearly, Eq. ([9]) indicates that the free energy of the system either increases or decreases as 
the magnetization of the free FDL (M2) rotates from the parallel (9 = 0) to the antiparallel 
{9 = tt). Since a lower free energy is favored, it means that M 2 is inclined to take one of 
the orientations depending on the sign of AQ n (7c). This preference on the magnetization 



orientation can also be described in terms of the magnetic energy of the FM layer when 
it is subject to a magnetic field: i.e., VLm{0) = — Heb-M-2 = —H e \,M. 2 cos6 ) where the x 
axis is chosen as the reference direction for the magnetic field and M. 2 — M 2 A tp is the 
total magnetic moment of the top FDL with A Q and tp its area and thickness. Through 
comparison with Eq. one can readily deduce the strength of the effective magnetic field 
that is essentially the exchange bias field mediated by graphene, 

(w) _ AO n (7r) 

Heb ~ 2M 2 ' 1 ' 

As in the case of conventional exchange bias, the strength of H$ is inversely proportional 
to the thickness tp [l| . However, one distinguishing feature of the exchange bias field given 
in Eq. (flOl) is its dependence on the electronic properties of the graphene layer, particularly 
the position of the electro-chemical potential \i that can be readily modulated by the gate 
bias (V g i, V g2 ; see Fig. 1). This also leads to qualitatively different characteristics for MLG 
and BLG that the calculation of AQ n (7r) highlights. 

Firstly, the signs of and H^' are different at least in the range < O.371. While 

MLG tends to establish M 2 parallel to Mi (H^ > 0), BLG favors the antiparallel alignment 

(2) 

{H eb < 0) as pointed out in our discussion earlier. Secondly, a shift of /i from the graphene 

charge neutrality point (e = 0) affects the strengths of the exchange bias fields in the 

opposite directions. Namely, the magnitude of gradually increases with \fi\, whereas 
(2) 

that of H eb is at the maximum at fi = and decreases to zero. A similar pattern is also 

observed in the temperature dependence. As T goes up, becomes stronger and 

weaker when /1 is at or around the graphene charge neutrality point. The properties of 

can be readily explained by considering the concentration variation of the conduction 

electrons and valence holes that can be spin polarized and therefore establish an effective 

field affecting equally the bottom and top FM layers. Hence, the parallel alignment of two 

ferromagnets are energetically favored and the magnitude of the indirect exchange bias field 

increases with the mediating carrier density (i.e., when MLG is used. On the other 

hand, the effect of H eb can be thought in terms of the band structure modification such as 

the size of the gap. As a larger gap pushes the energies of the occupied states lower, the 

(2) 

thermodynamic potential decreases likewise. This prefers a larger ang le 6 with H% < 
(i.e., antiparallel). Increasing and T diminishes the effect and, thus, the magnitude of 



H. 



(2) 



The aforementioned characteristics can be captured by an approximate expression in the 
limit of small G, 

A{l n (ii)~NG 2 f n (fi,T), (11) 

where N is the number of graphene primitive cells at the interface and the factor T) 
provides the specific dependence on \x and T for MLG (n = 1) and BLG (n = 2). The 
dependence of the exchange bias field [Eqs. ffTUj) and ffTTj) ] on the square of the effective field 
G is not surprising as the indirect interaction between the bottom (Mx) and the top (M 2 ) 
ferromagnets involves two interfaces with graphene. Figure 3 shows / n (/i, T) vs. evaluated 
at three different temperatures. As expected, the temperature factor has a considerable 
influence around fi = but its role diminishes significantly for \fi\ > 0.27! (7! = 0.4 eV). 
More crucial is the possibility of controlling f n (/j,T) (thus, Hj£') over a wide range even 
at room temperature. As evident from the figure, the shift of \x ~ ±0.157i can change 
the strength of by about a factor of two for both cases mediated by MLG and BLG. 
Based on these results, a rough estimate can be made for the magnitude of Assuming 
\f n \ ~ 1 meV (a typical number from Fig. 3) 20], a high temperature FM dielectric such as 
yttrium iron garnet, and G = 0.1 for the graphene interface (71 G = 40 meV) 9J, one can 
find tpM x ~ 1100 Oe, where tFM is the FM layer thickness in nm (e.g., H { 2 ] « 1100 Oe 
at tpM — 1 nm). Electrical control of the exchange bias in the estimated range can be of 
practical importance in a number of spintronic applications. 

In summary, our theoretical analysis shows that graphene can mediate the indirect inter- 
action of magnetic layers resulting in an effective exchange bias. Through the dependence 
on the graphene electro-chemical potential, it is also clearly identified that the effective 
exchange bias can be modulated electrically over a wide range even at room temperature. 
The numerical estimation indicates the potential significance of the proposed phenomenon 
in practical applications. 

This work was supported in part by the US Army Research Office and the FCRP Center 
on Functional Engineered Nano Architectonics (FENA). 
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FIG. 1: Schematic illustration of graphene sandwiched between two FM dielectric layers FM\ and 
FM2 of magnetization Mi and M2. While Mi is pinned by an AFM layer along the x direction, 
M2 can rotate on the x-y plane with 9 specifying the angle between them. The structure can be 
placed between metallic contacts providing electrical control with the gate voltages V g \ and V g 2- 
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FIG. 2: Energy band diagram of (a,b) mono- and (c,d) bilayer graphene when sandwiched between 
two FM materials as shown in Fig. 1. For (a) and (c), the angle 6 between Mi and M2 is zero 
(i.e., parallel alignment) and no gap exists between the bands. In the cases of (b) and (d), = tt 
(antiparallel) and the bands are doubly degenerate. The shaded regions schematically represent 
the occupied valence bands with fi = (dashed line) at T = 0. The effective field G = 0.1 is 
assumed for the calculation and the momentum is measured from the K or K 1 valley centrum. 
71 = 0.4 eV and vf = 10 8 cm/s. 
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FIG. 3: Factor f n (fi,T) [~ AU n (7r)/NG 2 ] vs. evaluated at three different temperatures (n = 1 
for MLG; n = 2 for BLG). This factor essentially determines the dependence of exchange bias field 
on fj, and T. 
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